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ABSTRACT 


Abstract of Report Presented to the Graduate School 
of the University of Florida in Partial Fulfillment of the 
Requirements for the Degree of Master of Science 

EVALUATION OF SEDIMENT TRAP EFFICIENCY IN AN ESTUARINE ENVIRONMENT 

By 

Daniel Mark Stoddard 
December 2001 

Chairman: Ashish J. Mehta 

Major Department: Civil and Coastal Engineering 

Trench-traps are utilized where sediment containment is a concern. In this study, trapping 
efficiency is key concern. A 60 m(L) x 300 m(W) x 2 m(D) trap was incorporated into the Cedar 
River, near the confluence with the Ortega River. A second trap of same dimensions was also 
incorporated 420 m upstream. Trap efficiency was calculated as a sediment removal ratio, or the 
percentage by which influent sediment load to the trap is reduced in the effluent load from the 
trap. Trap efficiency was carried out for varying Cedar River discharges. A specific discharge 
(16.4 m 3 /s) was found to yield the maximum removal. At discharges above and below this 
discharge, the removal ratio decreases. This is attributed to the increase in tidal influence at 
lower discharges and velocities too large to allow settling at higher discharges. 

Future work includes developing a monitoring scheme to determine actual sedimentation 
rates in a test trap at the chosen location. 
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CHAPTER 1 
INTRODUCTION 

1.1 Problem Statement 

Sediment shoaling in estuarine environments can create significant problems such as 
decreased discharges, degradation of water quality, and concentration of contaminants and 
organic compounds. One commonly employed solution to reduce sedimentation is the 
implementation of a trap scheme by creating a trench along the submerged bottom. To create a 
trench-trap, the depth at the chosen location is increased by dredging. In this study, a sediment 
trap is defined as an area of the submerged bottom deepened to a depth greater than the 
surrounding bottom, in order to reduce flow velocity. The lower velocity should allow sediment 
to deposit in the trap rather than move past and deposit elsewhere. This in turn allows for 
maintenance dredging to be performed at a specific location (the trap) rather than over a broad 
submerged area. The increased depth results in a decreased flow velocity, thereby allowing 
incoming sediment to settle in the trap itself, instead of being carried further downstream. The 
sediment can then be removed from the trap, rather than dredging the otherwise distributed 
deposit from a broader area. By holding the trap depth and location constant and varying the 
discharge of the river system, the efficiency of a trap can be assessed for different flow 
discharges. For present purposes, efficiency will be determined by the sediment removal ratio, 
which is the percentage by which the effluent sediment load (leaving the trap) is reduced with 
respect to influent load entering the trap (Ganju, 2001). By creating efficient traps much of the 
detrimental effects of excess sediment and unwanted pollutants entering the system can be 
curtailed. 
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1.2 Role of Florida Sediment 


In Florida’s biologically highly active estuarine and lacustrine environments, the fraction 
of fine-grained sediment that is organic is often on the order of 20-60% by weight and sometimes 
as high as 90-95%. There are three main sources of this organic matter. Terrestrial systems tend 
to be abundant in carbon (C), and the biomass produced by woodland and grassland is on the 
order of 50g C/m 2 (Mehta et al., 1997). Much of this material is degraded by the soil but some of 
it is washed away and introduced into fresh water and marine environments. The composition of 
this material is mainly cellulose which is non-degradable by water itself and the existing soil is 
less efficient in degrading the organic material making its breakdown very slow. Although 
aquatic plants breakdown more easily, they also contribute to the input of organic matter. The 
third source of organic matter is provided by phytoplankton, which usually has a biomass of 
1.5C/m 2 with 5-6 crops/year for the Florida region. Trefry et al. (1992) state that the coastal 
waterways in Florida are stressed by inputs of fine-grained organic-rich sediments from riverine 
systems. Besides the alterations of the benthic community that this input causes, there are 
indirect problems associated with organic sediment such as sorption of contaminants like Cd, Cu, 
Hg, Pb, Zn and PCB’s. In the Cedar River, PCB’s in sediment have been documented to be up to 
0.023 ppm (Campbell et al., 1993) and detectable amounts (up to 0.055 ppm) are also found in 
every species of fish collected from the area. 

Most of Florida estuaries are microtidal, hence another important feature of the region is 
the occurrence of episodic events such as heavy rainfall and storms which act as natural dredging 
mechanisms due to the strong currents generated (Marvan et al., 2001). 

The area of study is the Cedar River system located in Northeast Florida. Trapping 
contaminants in the Cedar River system is important due to the elevated concentration of PCB 
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(polychlorinated biphenyl) contaminants in the water system due to leeching of sediments and 
runoff from a fire at a chemical company in January 1984. The site was located approximately 
0.35 miles east of the Cedar River near the headwaters north of Interstate 1-10 and also adjacent 
to municipal storm drains and drainage ditches. The fire destroyed several tanks storing high 
concentrations (4,425 ppm) of PCB laden oils and other materials. It is believed that a 
combination of the damage to the storage tanks and the fire fighting effort created a vehicle for 
the PCB contaminant to enter the Cedar River basin. The surrounding groundwater and soil was 
sampled extensively in 1989 and the concentrations were still significantly above the regulated 
amount of 50 ppm. 

The filtering role of estuaries makes them crucial transitional areas trapping significant 
quantities of particulate and dissolved matter through a wide variety of physical and 
biogeochemical processes. Cohesive sediments play an important role in these processes. Unlike 
sand, well characterized by its grain size distribution, cohesive sediments are complex mixtures 
of different clay minerals, mainly organic matter, and a small percentage of sand and silt. 
Hydrodynamic action is the most important mechanism involved in sediment transport. It 
advects the suspended sediments, provides the force needed to erode the bed and, through 
turbulence, plays a major role in the flocculation of cohesive sediments. Relatively large 
velocities generally occur in tidal estuaries. Because the hydrodynamic processes involved in 
sediment transport are mainly non-linear, the sediments are very mobile in these estuaries. They 
are eroded and transported upwards during flood, deposited during slack water, eroded again and 
transported downwards during ebb and redeposited during next slack water, to restart their 
movement in the forthcoming tidal cycle. Cohesive and non-cohesive sediments are different 
from each other in two major aspects: flocculation and consolidation of deposited material with 
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compaction of the sediments. Floes are formed by joining individual particles and can strongly 
modify the settling velocity of particulate matter. After bottom deposition, the water content is 
still a significant part of the bed material. The expulsion of this water is part of the sediment 
consolidation process. The small pore dimensions imply long times for sediment deposition, 
which creates conditions for fluid-mud formation in environments with very high availability of 
sediments (Cancino and Neves, 1999). 

Fine-grained cohesive sediments are important in two types of engineering problems. The 
first relates to the sedimentation of harbors and channels and to dredging and navigation, and the 
second to the mixing and dispersion of contaminants. The properties of muddy sediments are 
significantly affected by chemical and biological factors. As a result of their cohesive nature, 
mud particles absorb pollutants, especially heavy metals and pesticides. As a result, 
understanding pollutant dispersion depends on an understanding of particle transport. The 
ubiquitous bacteria and other organisms secrete films that act as a very effective glue in 
enhancing the resistance of the bed to erosion. However, it is noteworthy that, in general, 
bioturbation acts both to increase cohesiveness and also to loosen beds and resuspend sediment 
(Mehta and Dyer, 1990). 

1.3 Objective and Tasks 

The main objective is to determine how sediment traps respond to variable discharges 
(concentration and velocity) in the Cedar River estuarine system by analyzing the trapping 
efficiency. 

Several tasks must be undertaken to determine the efficiency of this selected trapping 
scheme. These include: 
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1) Utilization of data collected from the field and from the existing literature to characterize 
the nature of the flow. This includes the tidal elevations and current velocities in the river 
system, and streamflow data for major tributaries from the literature. 

2) Analysis of the data collected to characterize the nature of the sediment to determine the 
historical suspended sediment concentration data. 

3) Modeling the flow field via a hydrodynamic model, in order to determine the velocities 
as well as the water surface elevations. 

4) Reevaluation and recalculation of rating curve results from previous analysis. This new 
curve will be used to calculate concentrations associated with the varying discharges. 

5) Utilization of a sediment transport model to determine suspended sediment 
concentrations. This model will incorporate the sediment characteristics determined from 
the sediment analysis. 

6) Two trap locations will be evaluated in the calibrated flow model, and the output from 
that model will be applied to the calibrated sediment transport model. The influent and 
effluent sediment loads through the trap will be recorded in order to quantify trap 
efficiency of each trap. 

1.4 Report Outline 

The following sections of the report will describe how the trap efficiency will be 
evaluated and the modeling efforts required. Next, a basic description of the Cedar and Ortega 
River system will be provided. The report will continue with development of the sediment rating 
curves, the trap design selection and the efficiency analysis of the selected trap. The summary, 
conclusions and recommendations will complete the report followed by a bibliography. 
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CHAPTER 2 

METHOD OF ANALYSIS 

2.1 Trap Efficiency 

The modeling of trap efficiency requires the use of a flow and sediment transport model. 
A flow model will provide water velocities and surface elevations, and these solutions will be 
applied to a sediment transport model. The sediment transport model will predict erosion, 
deposition, and suspended sediment concentrations in the presence of the trap. 

2.2 Flow Modeling 

2.2.1 Governing Equations 

The Navier-Stokes equations govern the free surface flows of constant density and 
incompressible fluids (Pnueli and Gutfinger, 1992). Applying the hydrostatic pressure 
distribution assumption yields three-dimensional model equations, and these can be vertically 
integrated to produce the following two-dimensional shallow water equations (Casulli, 1990): 
x-momentum: 


d(HU) djHUU) djHUV) 
dt dx dy 


„ drj d 
= -gH —- + — 
dx dx 
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continuity: 


( 2 . 2 ) 


^ + W + « = 0 (2 .3) 

dt dx dy 

where H is the water depth, U is the vertically-averaged horizontal x-direction velocity, V is the 
vertically-averaged horizontal y-direction velocity, t is time, g is the acceleration due to gravity, 
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r| is the water surface elevation measured from the undisturbed water surface, v is the eddy 
viscosity, and y is the bottom friction dependent coefficient defined as 



(2.4) 


where Ub and Vb are the horizontal x and y bottom velocity components respectively, and C z is 
the Chezy discharge coefficient, which is related to Manning’s n by 



(2.5) 


Solving this system of three partial differential equations (Equations 2.1, 2.2, and 2.3) for the 
three unknowns (U, V, q) can be accomplished via a numerical method. The numerical algorithm 
used is based on the method developed by Casulli (1990). First, a characteristic analysis is 
performed on Equations 2.1-2.3, in order to determine which terms must be discretized 
implicitly, such as the water surface elevation (Eqs. 2.1, 2.2), and the velocity divergence 
(Equation 2.3). The advective terms are discretized explicitly using an upwind scheme, which is 
unconditionally stable when a Eulerian-Lagrangian method is used to discretize the terms. This 
method requires the solution of a 5-diagonal matrix at every time step. It is used in conjunction 
with an alternating-direction implicit (ADI) routine, which results in two simpler, linear tri¬ 
diagonal matrices (Casulli, 1990). 

2.2.2 Model Operation 

The 2-D vertically averaged hydrodynamic model reported by Marvan (2001) which was 
developed by Casulli (1990) used in this study is operated using the MATLAB computational 
computer application. The use of MATLAB allows for the generation of the necessary graphics 
and data output in a simple fashion, though the computational effort is intensive, due to the 
necessity of large matrices. Rectangular grids with square elements are used, with numeric 
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“ones” indicating the body of water, and “zeros” representing land boundaries developed for 
input into the computer model. A similar grid is required for the input bathymetry, with the depth 
at mean high water entered into each element. 

2.2.3 Flow Boundary Conditions 

Flow boundaries are indicated by extending water cells to the grid edge. If freshwater 
inflow is desired, a permanent velocity can be imposed at the edge, corresponding to the desired 
flow condition. If a non-steady state inflow is desired, velocity as a function of time can be 
implemented. For a tidal flow boundary, a function specifying the water surface elevation at the 
boundary can be applied. If no velocity or elevation is specified at cells, which terminate at the 
grid edge, they become no-flow boundaries in the algorithm. 

2.2.4 Flow Model Input/Output Parameters 

The area and bathymetry grids described in Section 2.2.2 are required to specify the 
domain to be modeled. Other required inputs are the tidal forcing function at the seaward 
boundary, the calculation time step, the total simulation time, a file containing Manning’s n 
coefficient values for each cell, and velocities at the tributary flow boundaries. The output is 
three matrices consisting of the water surface elevations, x-direction velocities, and y-direction 
velocities, for every time step in the simulation. 

2.3 Modeling Sediment Removal 
2.3.1 Flow Model Input/Output Parameters 

The accurate prediction of suspended sediment transport in estuaries is important for 
activities such as dredging, the accurate mapping of navigation channels and improved 
understanding of pollutant transport. The difficulty in formulating accurate suspended sediment 
transport formulas suitable for a range of input conditions arises because the sediment transport 
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forms a complex feedback system with the near-bed hydrodynamics and the bed topography. 
One important area of research related to the accurate prediction of suspended sediment transport 
is the formulation of the magnitude and shape of the temporally averaged suspended sediment 
concentration profile in a tidal channel. This subject is of practical importance, as the product of 
the temporally averaged suspended sediment concentration and horizontal velocity profiles form 
the dominant component of the horizontal suspended sediment flux in tidal and steady flow 
conditions (Rose and Thome, 2001). 

The sediment removal will use a 2D MATLAB based horizontal depth-averaged fine 
sediment transport model with the capability to manipulate erosion and deposition functions 
based on organic content. Initially the model hydrodynamics are determined by a finite 
difference semi-implicit algorithm developed by Casulli (1990). In this, the water surface 
elevation is obtained implicitly and the velocity is determined in an explicit fashion. 

2.3.2 Governing Advection-Diffusion Equation 

Advection-diffusion is calculated with Equation (2.6) using a finite-volume explicit 
method based on the quadratic upstream (QUICKEST) method of Leonard (1977) (Marvan et al., 


2001 ): 


dhC | djhuC) | d(hvC) 8 
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dx 
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Dyx dx +Dyy dy 
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where t is time, C is the depth-averaged suspended sediment concentration, u and v are the 
longitudinal and transversal depth averaged velocities, h is the water depth and S is a source-sink 
term. The dispersion coefficients D xx , D xy , and Dyy are treated as follows (Preston, 1985): 
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( 2 . 8 ) 
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where AT/ and K t are the dispersion coefficients in the longitudinal and transversal directions 
taken to be 13 and 1.2 respectively, C z is the Chey coefficient and g is the acceleration due to 
gravity. The source-sink term in Equation (2.6) accounts for the erosion and deposition in the 
following way: 


S = Q e +Q d (2.10) 

where Q e is the erosion flux at every time step and Q d is the corresponding deposition flux. 

2.3.3 Deposition Flux 

The deposition flux is expressed according to Krone (1962) as: 

Q d =~pW s C (2.11) 


where W s is the sediment settling velocity and p is the probability for deposition defined as 

r T ' 

i-^ 


p = 


v 


( 2 . 12 ) 


l d J 


where r* the bed shear stress and r d a critical shear stress for deposition. In this analysis, the 
parameter r d is set to a value above the highest shear stress found in the modeled domain, thus 
allowing deposition to occur at all times as long as suspended sediment is present. 

Generally, the settling velocity of fine sediment is dependent on concentration. As a 
result, the settling velocity differs depending on three identifiable regimes: free settling, 
flocculation settling, and hindered settling. In the free settling range, relatively low 
concentrations permit the individual floes to settle without interference from other floes. The 
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settling velocity in this range is a function of the drag coefficient and the submerged weight of 
the floe. As concentration increases, the collision frequency of floes increases, resulting in the 
formation of larger floes. These floes are able to settle quicker due to their increased mass, and 
characterize the flocculation settling range. Eventually, the concentration in the water column 
reaches a point where a floe is unable to settle quickly due to significant interference from other 
floes, and the limited pore space for the fluid. This interference reaches a maximum when the 
water column resembles a bed of mud with negligible settling (Mehta, 1994). Hwang (1989) 
formulated a fit of the flocculation and hindered settling ranges, relating settling velocity to 
concentration as follows: 


aC 


(c 2 +»T 


(2.13) 


where a, b, m, and nw are empirical constants. For the free settling range, at concentrations 
(C<0.25 kg/m 3 ) a constant settling velocity (W s f) is provided. Laboratory tests performed in a 
settling column are required to determine the site-specific constants. The settling velocity of the 
aggregates is a function of concentration and of shear stress, because the aggregation of particles 
depends on the number concentration of particles in the flow, but their ultimate size is limited by 
turbulent shearing (Mehta and Dyer, 1990). 

2.3.4 Suspended Sediment Boundary Conditions 

The boundary conditions at the tributary connections can be expressed as steady-state 
concentrations, or sediment rating functions can be applied if unsteady tributary flows are 
desired. This also holds true at the tidal entrance, where incoming concentrations can be 
specified, varying with tidal stage and/or current velocity. 
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2.3.5 Erosion Flux 


In contrast to tj, Q e cannot be treated in this way because erosion depends on the shear 
strength of the soil and is therefore considered in the following form: 

Qe= £ A^~f s ) (2-14) 

where £n is the erosion rate constant and r s is the bed shear strength. The bed shear stress is 
computed as 


(#+# 


where pw is the density of water. The shear strength of the bed is calculated via Mehta and 
Parchure (2001): 

< 2 - 16 ) 

where tf is the solids volume fraction (pr/p s ), Pd is the dry density, p s is the grain density, <f> s is 
the limiting solids volume fraction value of ^at which zj = 0 and a and (3 S are sediment-specific 
empirical coefficients. 

2.5 Sedimentation, Sediment Trap, and Trap Efficiency 


2.5.1 Sedimentation 

Sedimentation at any point in the estuary can be calculated from the deposit thickness Ah 
given by 


& O. A t 
M Pd 


(2.17) 


where t s is the total simulation time, At is the time step, i is the time step index, Qdi is the 
deposition flux, and pa is the deposit dry density. The sedimentation rate is then calculated by 


12 



(2.18) 



where Sr is the sedimentation rate. 

It is well known that in rapidly moving waters, fine particles are not deposited. 
Consequently, consideration must be made that sedimentation increases as flow velocity 
decreases. Since there is a close relationship between flow velocity and bottom shear, the 
following hypothesis has been proposed: shear stress in moving waters is an important 
controlling (reducing) factor on sedimentation. In aquatic systems, there are cohesive and 
noncohesive particles. The theoretical basis of cohesive material sedimentation in moving waters 
was derived from the fundamental flume experiments of Partheniades (1965, 1972). These 
experiments showed that sedimentation, as well as resuspension, depend on bottom shear stress. 
There are threshold values, one for sedimentation and a much higher value for resuspension. It 
was found that the degree of deposition (the proportion of the initially resuspended material, 
which settles) as well as the rate of deposition is controlled by the bottom shear. Partheniades 
(1972) wrote: “We may distinguish two groups of floes: those with sufficient high strength to 
resist the flow induced disruptive shear stresses, which are highest near the bed, and those with 
insufficient strength. The first will be able to reach the bed, will develop several bonds with it 
and will become a part of it; the remaining floes will be disrupted and reentrained” (Kozerski and 
Leuschner, 1999). 

2.5.2 Definition of Trap 

In this study, a sediment trap is defined as an area of the submerged bottom deepened to a 
depth greater than the surrounding bottom, in order to reduce flow velocity. The lower velocity 
should allow sediment to deposit in the trap rather than move past and deposit elsewhere. This in 
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turn allows for maintenance dredging to be performed at a specific location (the trap) rather than 
over a broad submerged area. 

2.5.3 Definition of Trap Efficiency 

Trap efficiency is defined as the percent by which effluent suspended sediment load is 
reduced with respect to the influent suspended sediment load (removal ratio). In a tidal situation, 
the seaward edge of the trap will be the influent side during flood tide, and the effluent side 
during ebb tide, and vice versa for the landward edge. 

2.5.4 Calculation of Trap Efficiency 

At each time-step in the sediment transport simulation the concentration, the velocity, and 
the water surface elevation will be calculated in each cell. The cells that border the trap and are 
flow-normal are also of interest. Sediment loads can be calculated for these border cells as 
follows: 

q - UCHAx (2.19) 

where q is the sediment load, U is the flow velocity, and Ax is the cell width. The sediment load 
on each side of the trap will be used to compute the sediment removal ratio as follows: 

R=‘LZQjl (2.20) 

<h 

where R is the removal ratio, q t is the influent sediment load, and q e is the effluent sediment load. 
The removal ratio will be averaged over a tidal cycle, using the removal ratio values from each 
time-step. 

2.5.5 Calculation of Trap Efficiency as a Function of Discharge 

Simulations will be run for different discharges of the Cedar River. The average removal 
ratio for each trap will be compared for a given flow discharge. The normal flow will be used as 
the benchmark by which the efficiency of the trap will be assessed under the other discharge 
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cases. The removal ratio as a function of discharge will be plotted to determine what effect 
discharge has on trap efficiency. 

CHAPTER 3 

CEDAR, ORTEGA, AND ST. JOHNS RIVER SYSTEM 
3.1 History and Description of the System 

Before European involvement in North America, the Timucuan Indians called the St. 
Johns River Welaka, or river of lakes. In the early 1500s, Spanish seamen called the river Rio de 
Corrientes or River of Currents. In 1562, almost 50 years before the settlement in Jamestown, the 
French established Fort Caroline on a high bluff overlooking a river they called Riviere de Mai 
(River of May) because they arrived there on May 1. In 1565, Spanish soldiers marched north 
from St. Augustine, captured Fort Caroline and slaughtered the French. The Spanish renamed the 
river San Mateo to honor the saint whose feast followed the day they captured the river. Later, 
the river was renamed Rio de San Juan after a mission near its mouth named San Juan del 
Puerto. The English translation of the name Rio de San Juan, St. Johns River, lasted through 
English, Confederate and American possession of the river and remains today. Soon after 
England acquired Florida in 1763, King George III sent botanist John Bartram to explore Florida. 
His son, William Bartram, stayed in Florida and published his book Travels in 1791. It describes 
his exploration of the river as far south as Lake Harney. In the 1800s, steamboats made the St. 
Johns River a popular winter destination for northerners. By the 1860s, several steamers were 
making weekly round trips from Charleston and Savannah to Jacksonville and Palatka, and other 
settlements. In the 1900s, miles of floodplain were drained to make room for indigo, sugar cane, 
citrus and other profitable crops. Encroachment through draining of the headwater marshes at the 
river’s southern end was neither planned nor controlled. More than 70 percent of the marsh was 
claimed for agricultural and urban uses. In 1954, Congress authorized flood-control works in the 
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southern part of the St. Johns River. To store water and to move floodwaters, large reservoirs and 
canals were designed by the U.S. Army Corps of Engineers. The Corps’ project was halted in the 
1970s. In 1974, the project was deemed unacceptable for environmental reasons. In 1980, a 
redesigned project by the St. Johns River Water Management District favored restoring wetlands 
to hold and release floodwaters and managing water levels to simulate natural marsh conditions. 
Since the project began, the District has restored more than 610 km 2 of original marsh, an area 
about the size of Delaware. 

The St. Johns River is an ancient intracoastal lagoon system. As sea levels dropped, 
barrier islands became an obstacle that prevented water from flowing east to the ocean. The 
water collected in the flat valley and slowly meandered northward, forming the St. Johns River. 
The St. Johns River is the longest river in Florida at 500 km in length. The width of the river 
varies between a flat marsh at its headwaters and averages about two miles in width between 
Palatka and Jacksonville. In central Florida, the St. Johns River widens to form large lakes. 
Additionally, it is one of the few rivers in the United States that flows north. The total drop of the 
river from its source in marshes south of Melbourne to its mouth in the Atlantic Ocean near 
Jacksonville is less than 9 m, or about 1.6 cm per kilometer, making it one of the “laziest” rivers 
in the world. Because the river flows slowly, it is difficult for the river current to flush pollutants. 
Major pollution sources include discharges from wastewater treatment plants and stormwater 
runoff from urban and agricultural areas. This runoff carries pesticides, fertilizers and other 
pollutants into canals, ditches and streams that lead to the river with much of the river pollution 
is concentrated around urban areas. Salt water enters the river at its mouth in Jacksonville. In 
periods of low fresh water flow, tides may cause a reverse flow as far south as Lake Monroe, 260 
km upstream from the river’s mouth. 
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The Ortega River basin is located west of the St. John’s River in south-central Duval 
county in northern Florida and is an important tributary of the St. Johns River (Figure 3.2). The 
Ortega River is the main tributary of the system, discharging approximately half of the total 
system’s volume to the St. Johns River. The Cedar River is the second most important tributary 
and there are three other secondary tributaries of the system (Fishing Creek, Butcher Pen Creek 
and Williamson Creek). The upstream portion of the Ortega River is known as McGirt’s Creek. 
The creek lies within the Duval uplands physiographic province and flows generally north to 
south. The Ortega River continues this course until it reaches the Eastern Valley physiographic 
province, where the river gradually turns 180 degrees to a north-northeasterly course before 
reaching the St. John’s River north of the Jacksonville Naval Air Station. 







Figure 3.2: Cedar/Ortega River system and tributaries (from Mapquest.com). 

The Cedar River, the largest tributary of the Ortega River (Figure 3.2), is actually a major 
system itself. From its head waters north of Interstate 10 and west if Interstate 295, this river 
flows southeast to its confluence with the Ortega River. Major tributaries to the Cedar River are 
Willis Branch, Williamson Creek, Butcher Pen Creek, and Fishing Creek. The tidal interface for 
the Ortega River is at Collins Rd., while the tidal interface for the Cedar River is near Lane Ave. 
(These two and other road locations are not highlighted in any drawings herin; they are found in 
road maps of the Jacksonville area.) 

In the early 1990’s, approximately one-third of the Ortega/Cedar River basin was residential, 
with commercial/industrial and vacant land comprising the other major land uses. Since then, 
vacant land has decreased significantly (Mehta et al., 2000). The average annual rainfall in the 
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Ortega/Cedar basin is approximately 132 cm and the major portion of it falls between June and 
September (Campbell et al., 1993). Water depths in the Ortega/Cedar basin study area range 
between 1 and 7 m, with the range in the Cedar River between 0.65 m and 4.3 m. At the mouth 
of the Ortega River with St. Johns River, the semidiurnal tide ranges from 0.14 m (neap tide) to 
0.28 m (spring tide) having a mean of 0.19 m (this was the value used in the hydrodynamic 
model). The bottom and suspended sediment is mostly a mixture of clay, silt and organic matter. 
Typical suspended sediment concentration is approximately 15 mg/1; however, during storm 
runoff events it rises to as much as 105 mg/1. The mean organic content was found to be 28%. 
Previous samples from Mehta et al. (2000) show similar results for the sampled area having 
values between 22 and 36%. Measurements were also obtained from the St. Johns River Water 
Management District (SJRWMD), which showed less organic content ranging from 8-22% 
within the study area (Marvan et al., 2000). For this analysis the mean organic content value will 
be used for calculation purposes. 


19 




CHAPTER 4 

ASSESSMENT OF TRAP EFFICIENCY 


4.1 Sediment Rating Relations 
4.1.1 Sediment Rating Curve Definition 

In many lowland rivers a major part of the sediment is transported in suspension. As the 
finest fraction of the suspended sediment load often is a non-capacity load it cannot be predicted 
using stream power related sediment transport models. Instead, empirical relations such as 
sediment rating curves often are applied. A sediment rating curve describes the average relation 
between discharge and suspended sediment concentration for a certain location. The most 
commonly used sediment rating curve is a power function (e.g. Walling, 1974, 1978): 

C = aQ b (4.1) 

where C is suspended sediment concentration (kg/m ), Q is water discharge (m /s), and a and b 
are regression coefficients. Equation 4.1 covers both the effect of increased stream power at 
higher discharge and the extent to which new sources of sediment become available in weather 
conditions that cause high concentrations are related to the statistical method used to fit the 
sediment-rating curve and to the scatter about the discharge. Despite its general use several 
problems are recognized that regard the accuracy of the fitted curve as well as the physical 
meaning of its regression coefficients. Statistical inaccuracies related to the fitting procedure are 
discussed by Ferguson (1986, 1987), Jansson (1985), Singh and Durgunoglu (1989), and Cohn et 
al. (1992). They concluded that the sediment load of a river is likely to be underestimated when 
concentrations are estimated from water discharge using least squares regression of log- 
transformed variables. Scatter, among other things, caused by variations in sediment supply due 
to, for instance, seasonal effects, antecedent conditions in the river basin, and differences in 
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sediment availability at the beginning or the ending of a flood. This is not accounted for by the 
rating curve. As a sediment rating curve can be considered a ‘black box’ type of model, the 
coefficients a and 6 in equation 4.1 have no definite physical meaning. Nevertheless, some 
physical interpretation is often ascribed to them. Peters-Kiimmerly (1973) and Morgan (1995) 
state that the a-coefficient represents an index of soil erosion severity. High a-values indicate 
intensively weathered materials, which can easily be eroded and transported. According to 
Peters-Kiimmerly (1973), the ^-coefficient represents the erosive power of the river, with large 
values being indicative for rivers where a small increase in discharge results in a strong increase 
in erosive power and sediment transport capacity of the river. Others state that the ^-coefficient 
indicates the extent to which new sediment sources become available when discharge increases. 
Several authors compare the values of the ^-coefficient obtained for different rivers to discuss 
differences in sediment transport characteristics in the different basins (Peters-Kiimmerly, 1973; 
Walling, 1974; Sarma, 1986; Morgan, 1995; Kern, 1997). 

As discussed previously, the values of the regression coefficients of sediment rating 
curves are assumed to depend on the severity of erosion, or the availability of sediment in a 
certain area, the power of the river to erode and transport the available material, and on the 
extent to which new sediment sources become available in weather conditions that cause high 
discharge. According to Walling (1974) 6-values are also affected by the grain size distribution 
of the material available for transport, i.e. in streams characterized by sand sized sediments the 
power of the stream to transport sediment will be more important than in streams that mainly 
transport silt and clay. This will result in high 6-values. However, as the a- and 6-coefficients of 
sediment rating curves are inversely correlated (Rannie, 1978; Thomas, 1988) it seems more 
appropriate to use the steepness of the rating curve, which is a combination of the a- and 6- 
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values, as a measure of soil erodibility and erosivity of the river. Steep rating curves, i.e. low a- 
and high 6-values, should thus be characteristic for river sections with little sediment transport 
taking place at low discharge. An increase in discharge results in a large increment of suspended 
sediment concentrations, indicating that either the power of the river to erode material during 
high discharge periods is high, or that important sediment sources become available when the 
water level rises. 

Flat rating curves should be characteristic for river sections with intensively weathered 
materials or loose sedimentary deposits, which can be transported at almost all discharges. When 
this line of reasoning is accepted for the Cedar River, the following interpretation can be 
assigned to the rating curves shown in Figure 4.3. This suggests a limited amount of fine 
sediment, which can be picked up from the bed at low discharge. Once a certain discharge 
threshold is exceeded, sediment supply to the river increases, and sediment can be picked up 
from the riverbed, resulting in a rapid increase in suspended sediment concentrations. This 
argument leads to the hypothesis that steeper rating curves are indicative of rivers, or river 
sections, where most of the sediment transport takes place at high discharge. Sediment transport 
rates estimated using a sediment-rating curve always differ somewhat from measured sediment 
loads. Hence, the sediment-rating curve produces only reasonable estimates of long-term total 
sediment loads. (Asselman, N.E.M., 2000). 

4.2 Determination of Rating Curve 

In the original analysis performed by Marvan (2001), the number of suspended sediment 
measurements were limited, hence, in order to set the boundary conditions for the tributaries a 
rating curve was developed by relating the peak and mean values of concentration C to the 
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corresponding peak and mean river discharge Q in the Ortega and Cedar rivers, obtaining the 
following relations (Fig. 4.1). 

C = 2xlO~ 5 0 223 (4.2) 

At the mouth of the Ortega river, this boundary condition was only applied when the flow 
was entering the system from St. Johns River. Since no concentration data were available from 
this site, by using the sediment characteristics and a zero-dimensional resuspension model 
(Mehta and Li, 1999) the following rating curve was developed for the Cedar river: 

C = 1.65x10~ 2 (9° 52 (4.3) 

At the time of Marvan’s analysis, only limited data were available. Further analysis was 
performed to validate the Cedar River rating curve utilizing additional data from the SJRWMD. 
The additional discharge data were included with original data and averaged by month and 
plotted against the concentration, the resulting data were band averaged to reduce scatter. A new 
rating curve was fit to the band averaged data as shown in Figure 4.1. Following is the 
recalculated rating curve for the Cedar River: 

C = 6xl<r 2 0 028 (4.4) 
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Revised Cedar River Rating Curve 

(data band averaged over time by month) 



Figure 4.1 Revised Cedar River sediment rating curve. 


Band Averaged Rating Curve & Marvan (2001) 
Rating Curve (Cedar River) 



Figure 4.2 Comparison between Marvan (2001) and new Cedar River sediment rating curves. 
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Cedar/Ortega Sediment Rating Curves 



Ortega Mouth 

New Cedar 
River 


Figure 4.3 Cedar/Ortega River sediment rating curves. 

Figure 4.2 shows a comparison between the original Cedar River rating curve and the 
revised curve. The new curve values were incorporated in to the sediment model. Figure 4.3 
shows the sediment rating curves for the Cedar/Ortega system that serve as sediment 
concentration boundary conditions for the sediment transport model. 

4.3 Trap Design Selection 
4.3.1 Factors and Considerations 

Some basic design factors were considered in sizing and locating the traps. The following 
are some basic factors as Parchure et al. (2000) indicates to consider when designing a sediment 
trap: 

a. Locate the trap at a place of maximum sediment transport. 

b. It should have navigational access for a dredge to get in and get out without difficulty. 

c. The depth and size of the trap should permit safe operation of a dredge. 
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d. The storage volume of the trap should permit adequate temporary storage of the 
sediment. 

e. Preferably, the trap should catch both fine and coarse sediment. 

f. The prevailing flow pattern should be approximately normal to the longer side of the 
trap. 

These factors were evaluated and applied in varying degrees when selecting the size and location 
of the test traps. 

4.3.2 Evaluation 

Based on flow data and previous hydrodymanic and sediment transport analyses of the 
Cedar River system by Marvan et al. (2000) and Mehta et al. (2000), Trap 1 is placed near the 
confluence of the Cedar and Ortega Rivers. Trap 2, conversely is located approximately 420 m 
upstream. Both traps were selected 60m (1 cell) wide by 300m (5 cells) long with a surface area 
of 18,000 m 2 and a volume of 36,000 m 3 . The traps are to have an initial dredged depth of 2 
meters (from original bed depth) of the river cross-section, which for Trap 1 the dredging depth 
is 3.8 m and the dredging depth for trap 2 is 3.2 m. These were considered sufficient to reduce 
the velocity in the canal to allow measurable sediment to settle. For example, with a flow of 3 
m 3 /s and regular tidal forcing in the Cedar River, the average velocity over trap 1 was found to 
be 0.13 m/s. Over the same location with no trap in place, the average velocity was 0.24 m/s, 
which results in a 49% reduction in velocity over the trap. The range of velocity reduction for 
trap 1 and trap 2 were from 44% to 51% and 43% to 53%, respectively with an overall average 
velocity reduction of 48%. 
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Figure 4.4 Bathymetry of Cedar/Ortega River as used in hydrodynamic/sediment transport 
models. 



Figure 4.5 Bathymetry of Cedar River as used in hydrodynamic/sediment transport models. 
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In order to correctly incorporate the trap into the computer model, the existing 
bathymetry file was updated with the location and depths of the traps at the selected cells. The 
output from the flow model was then input into the sediment transport model. Sediment removal 
ratio, as defined by Equation 2.20, was calculated from the influent and effluent sediment loads 
in units of kg/s (Equation 2.19) at the cells adjacent to the trap on the upstream and downstream 
edges of the trap. The input densities (dry, bulk, granular) required for the trap were determined 
from Figure 4.6, and the values are 157 kg/m 3 , 1099 kg/m 3 , 2188 kg/m 3 , respectively. 



♦ Granular density 

* Bulk density 
Dry density 

+ Loxahatchee River 


0 10 20 30 40 50 60 70 


Organic content (%) 

Figure 4.6 Variation of granular, bulk, and dry densities with organic content using data from 
three Florida locations and the Loxahatchee River (from Ganju, 2001). 
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Figure 4.7 Settling velocity vs. sediment concentration. 


Fishing Creek 


Cedar/Ortega 
Confluence Area 



Figure 4.8 Cedar River section of the computational grid. (Trap cells are shown in black.) 
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Influent and effluent sediment loads were calculated for each time step, and then the 


removal ratio averaged over one ebb tidal cycle as follows: 

M 


R 


i-i 


M 


(4.5) 


where R ave is the ebb tide averaged removal ratio, i is the index for each time step At, R; is the 
removal ratio from a single time step, and M is the total number of time steps over a complete 
tidal period, as follows: 


M = 


T 

2At 


(4.6) 


where T is the tidal period (12.42 h). Flood tidal data were not used to calculate the removal ratio 
because the effluent load (downstream edge) contained only the sediment which escaped the trap 
on the previous ebb tide. It was observed that for 63% flow (35.4 m 3 /s) was the only discharge 
that produced a constant flow toward the St. Johns River even while the tide was flooding 
although small flood values appeared at a discharge of 16.4 m /s. For the 35.4 m /s flow, the 
removal ratio was calculated for the entire tidal period. 

4.4 Trap Efficiency as a Function of Discharge 
4.4.1 Trap Performance 

Harmonic analysis was carried out by Marvan (2001) for tide at the Ortega River mouth 
at the St. Johns River in order to generate a one-year water level record. Water surface elevation 
and velocity data from San Juan Bridge (Fig. 3.2) provided by SJRWMD was used for 
calibration of the hydrodynamic model. River discharge was also available in Ortega River 
having a mean discharge of 1.4 m 3 /s and a maximum of 112 m 3 /s. By measuring the watersheds 
of the other main tributaries (Fishing Creek, Butcher Pen Creek, Williamson Creek and Cedar 
River), an estimate of the river discharge was made, yielding the rates given in Table 4.1. 
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Table 4.1. - Cedar/Ortega and tributary discharges in m 3 /s _ 

Tributary _ Normal conditions _ Storm runoff event 

Ortega River 8.50x10' 7.80x10 

Fishing Creek 1.60x10'* 1.46x10* 

Butcher Pen Creek 4.00xl0‘ 2 3.75x10° 

Williamson Creek 3.90xl0' 2 3.64x10° 

Cedar River _ 6.50x10'* _ 5.52x10* _ 

(From Marv&n et al., 2001) 

The hydrodynamic and sediment transport models used were previously calibrated as 
discussed in Marvan et al. (2000) and utilized for the trap analysis. Using the calibrated model, 
several discharges were used for evaluation. Table 4.2 provides the evaluation discharges and 
associated concentrations using the recalculated Cedar River rating curve. Due to model 
limitations at the time of evaluation, the 100% discharge of 55 m 3 /s was not evaluated. The 
organic content included in the model was 28%, which is the mean organic content of the Cedar 
River sediment. Removal ratios were calculated only during periods of ebb tide flow through the 
trap (Section 4.3), and plotted against Cedar River discharge (Figure 4.9) and velocity (Figure 
4.10). These simulations show that the removal ratio is maximum at a discharge of 
approximately 16.4 m 3 /s. At higher discharges the removal declines meaning that the velocity is 
too large to allow particles to settle in the trap and are subsequently transported past the trap. 
Conversely, at significantly lower discharges the same particles settle before arriving at the trap. 

To provide a performance comparison, the trap was moved 420 m upstream from its 
previous location. The dredge depth and surface area of the trap remained the same at 18,000 m . 
The original bathymetric grid was adjusted to reflect the new depth. The removal ratio for Trap 2 
was calculated for the same discharges. Table 4.2 compares the removal ratios for the discharges 
at each trap location. Trap 2 performed 28 % better at the peak removal ratio flow rate (16.4 
m 3 /s) and by an overall average of 56% over Trap 1. This reduced performance by Trap 1 can be 


31 



partly attributed to the increased tidal action near the confluence of the Cedar and Ortega Rivers 
where Trap 1 is located. Trap 2 performed more effectively due to more consistent flow direction 
and velocity since the location is well within the Cedar River. Trap 1, being closer to the Ortega 
River confluence area experienced more “mixing”, reducing trapping efficiency. 

Table 4.2 Removal ratio as a function of Cedar River discharge for Trap 1 and Trap 2 


Cedar River 
Discharge (m 3 /s) 

Trap 1 Removal Ratio 

Trap Removal Ratio 

0.65 

.09 

.14 

3.0 

.18 

.29 

16.4 

.27 

.35 

34.5 

.15 

.26 


Removal Ratio vs. Cedar River Discharge 



Figure 4.9 Removal ratio of trap 1 and 2 as a function of Cedar River discharge. 
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Removal Ratio vs. Discharge Velocity 



-•-Trap 1 
—A— Trap 2 


Figure 4.10 Removal ratio of trap 1 and 2 as a function of Cedar River velocity. 

4.4.2 Tidal Influence on Performance 

Similar to the results in Ganju (2000), the removal ratios for both traps increase until a 
critical discharge is reached and then decrease. For this system, each trap, the “break-even” 
discharge was 16.4 m 3 /s. Considering the graphical removal ratio solution from Mehta (1984) for 
a basin with two entrances, the Cedar River trapping schemes results in the relationship 

Roc! (4.7) 

U 

with settling velocity, trap length constant, the removal ratio is proportional to the discharge 
velocity. As verification, recalling the depositional flux is: 

Q d =W s CL (4.8) 

and the total sediment flux entering the trap is 

Qi = CUh (4.9) 

and taking the removal ratio as 
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(4.10) 


q, a Uh 

it can be shown to be proportional to U" 1 since W s , L, and h are essentially constant. 
Conceptually, as the velocity decreases, the removal ratio should increase. This theory is 
supported by Baker et al. (1988), but for both of these situations, the velocity was unidirectional. 
Baker et al. showed that for a sediment trap in a unidirectional flow the removal ratio increased 
as the velocity approached zero. For the chosen trap locations, as the discharge velocities 
decrease, tidal influences increase. Starting with the highest discharge (unidirectional flow) and 
working toward the minimum discharge (bi-directional flow), the removal ratio increases until a 
maximum value is reached for each trap. Traps 1 and 2 had maximum removal ratios of 0.27 and 
0.35, respectively. Continuing to reduce the discharge from the peak removal ratio at 16.4 m /s, 
the tidal influences begin to appear changing the flow from unidirectional to bi-directional. As 
the discharge approaches zero the tidal forcing effect increases and becomes maximum, which 
keeps the sediment in the vicinity of the trap in a semi-resonant pattern for a longer period of 
time before being pushed through the system as in the unidirectional or runoff induced flow case. 

Based on the results and observations thus far, it is believed that tidal influence is a 
contributing factor in removal ratio calculation and should be evaluated. In a unidirectional flow, 
the sediment is more likely to settle because external disturbances are reduced versus the 
directional velocity changes that take place in a bi-directional flow situation. By changing 
direction, some of the previously deposited sediment may become resuspended if consolidation 
has not occurred. If the discharge is minimal and tidal forcing is near maximum, the 
consolidation would be small and resuspension more likely. With this in mind, the removal ratio 
would decrease. Of course, actual sediment characteristics that would have to be evaluated for 
each system depends significantly on settling velocity. In this analysis, the settling velocity is 
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free settling and constant because the concentrations are small. Each estuarine system would 
have to be evaluated to determine the discharge value where tidal influence begins to impact 
removal ratio. A general relationship can be developed for both tidal and non-tidal influenced 
removal ratio portions. 

As tidal influence increases, the equivalent trap length theoretically increases due to the 
resonance in the system. This tidal equivalent trap length (L) can be illustrated by using the 
following expression: 


L_ 

K 



m > 1 


(4.11) 


where L 0 is the original trap length and U 0 is a characteristic velocity of the non-tidal portion of 
the removal ratio (most likely for the desired evaluation discharge). The value m is a scaling 
factor to account for the varying tidal influences as the discharge changes. Substituting this 
expression into Equation 2.30 and converting the velocities into discharges, the following 
expression results for the tidal-influenced removal ratio: 
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where X is the dimensionless tidal influence factor. Figure 5.2 shows the tidal/non-tidal 

3 

influenced removal ratio using Equations 2.30 and 2.32 for k=l, m=1.5, and Q 0 = 5 m /s.. 
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Figure 5.1 Single-Box model for illustration of tidal/non-tidal removal ratio as a function of 
deposition and erosion fluxes. 
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Figure 5.2 Tidal/Non-tidal removal ratio as a function of discharge 

The resulting trend from Equation 4.12 is very similar to the model results verifying that a tidal 
influence that reduces the removal ratio is present at lower discharges. To ensure a valid and 
accurate application, the appropriate scaling factors (m and k) need to be adjusted for each 
system and trap configuration. 
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CHAPTER 5 
CONCLUSIONS 

5.1 Summary 

The objective of this study was to determine the effect of discharge on trapping efficiency 
for a given trap design in different locations. The Cedar River was chosen as the location of the 
study due to the influx of organic rich fine sediment and contaminants from upstream sources 
desired to be kept in the traps rather than spread through the entire biologically sensitive 
estuarine system. Flow (hydrodynamic) and transport models were utilized to calculate the water 
levels, velocities, and sediment concentrations in the Cedar/Ortega River system. The models 
were calibrated using data previously collected from available sources and field investigations 
(see Marvan, 2001). Tidal elevations and currents were measured, and historical tributary flow 
data were obtained from St. Johns River Water Management District (SJRWMD in order to 
calibrate the flow model. 

Utilizing the calibrated model, a sediment trap in two locations was incorporated into 
each of the models to determine the trapping efficiency as a function of discharge. The results of 
these simulations and conclusions are discussed in the following section. 

5.2 Conclusions 

1) The simulations for trap efficiency as a function of Cedar River discharge demonstrate a 
specific discharge (16.4 m 3 /s) at which the sediment removal ratio is a maximum. Above this 
discharge, particles are moving fast enough to bypass the trap and below this discharge the 
particles deposit before arrival at the trap. A comparison trap was evaluated 420 m upstream. 
The trap performed 28% better while maintaining the maximum removal ratio at 16.4 m /s. 
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2) Comparing the trapping efficiency results against the expected relationship for removal 
ratio, the inconsistencies appear at lower velocities. According to the theory, as velocity 
decreases the trapping efficiency increases since the removal ratio is proportional to U '. This 
indicated another influence was present in the system reducing the efficiency at lower velocities. 
Also noticed as the discharge velocities decreased, the tidal influence became stronger. A 
relationship was developed and applied to account for the increase tidal influence at lower 
velocities Equation 4.12. 

5.3 Recommendations for Further Work 

The trapping efficiency calculation would be more accurately performed with a 3-D 
model to more effectively account for the mud suspension in high concentrations just above the 
bed. 

Prior to any full scale dredging, a test pit should be dredged to accurately determine 
sedimentation, flow velocity (to calculate discharge), pressure (to determine water elevation), 
and deposition thickness in the test pit in the Cedar River. A proposed trap is showed in Figure 
5.1. The system setup would have a central data logger to record output from the turbidity 
meters, current meter, and pressure sensor. To determine the equilibrium bed elevation and the 
pit bed level, divers would be needed. A similar system was deployed by Harley and Dean 
(1982) off the coast of Colombia. 
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1. Onset {Bourne, MA) 
Tattletale 6 data logger 


3. Falmouth Scientific (Cataumet. MA) 
2*D ACM current meter 


2. Seapoint (Kingston, NH) 4. Trans Metrics (Solon, OH) 



Figure 5.1 Possible layout of an experimental test pit from Ganju (2001). 

Vicente (1992) provides a method to calculate a constant sedimentation coefficient (K), 
which can be used to calculate bed elevation through time, as follows: 

= (5-0 

where H(t) is the bed elevation at any time, H 0 is a datum elevation, and H e is the equilibrium 
bed elevation (in absence of dredging). In order to determine K, test pit bottom elevations must 
be recorded over time, referenced to H e . Once the time/elevation data are recorded for different 
areas of the pit, K is determined. One benefit of this method is that only two monitoring visits are 
needed to determine K. Shoaling depth through time can then be estimated for the site, with the 
final shoaling depth approaching the equilibrium bed elevation (Ganju, 2001). 
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